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Abstract. The last part of SpS5 dealt with the circumstellar environment. Structures are indeed 
found around several types of massive stars, such as blue and red supergiants, as well as WRs 
and LB Vs. As shown in the last years, the potential of IR for their study is twofold: first, IR can 
help discover many previously unknown nebulae, leading to the identification of new massive 
stars as their progenitors; second, IR can help characterize the nebular features. Current and 
new IR facilities thus pave the way to a better understanding of the feedback from massive stars. 
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1. Introduction 

Circumstellar material holds clues about the mass-loss history of massive stars. Indeed, 
as the winds interact with the interstellar medium (wind-blown bubbles, bow shocks), 
they leave a characteristic signature that depends on the wind properties. Moreover, the 
material ejected during short eruptive phases is visible as nebulae around massive stars. 
The analysis of these features reveals which material was ejected and in which quantity. 
With the recent reduction in mass-loss rates, these episodes of enhanced mass-loss have 
gained more attention, as they seem more crucial than ever in the evolution of massive 
stars. 

Another reason to study the close environment of massive stars is to better understand 
the evolution of supernova remnants (SNRs). Indeed, the famous rings of SN1987A may 
only be understood if one considers the previous mass-loss episodes of the progenitor. 
Morphology is not the only SNR parameter which is affected, as the SNR dynamics in 
an homogeneous medium or in winds and circumstellat eject a is not identical. 

For its study, the IR provides several key diagnostics. Continuum emission in this range 
is provided by heated dust, which may have a range of temperatures depending of the 
framework (very close hot features, large, old, and cool bubbles). In addition, IR lines 
probe the many phases of the material: molecules (e.g. PAHs) for the neutral material, 
ionized metals for HII regions,... 

This summary of SpS5 - part III examines each case of circumstellar environment in 
turn, and concludes with the potential offered by current and future facilities. 

t FNRS Research Associate 
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Figure 1. Left: Ha emission (greyscale) of Vela X-1 with PACS 70/im emission contours shown 
on top. Right: Colour composite image of bow shock of BD+43°3654 (WISE 12/im in blue, 
PACS 70/im in green, and PACS 160/im in red). The direction of proper motion is indicated by 
the arrow in both cases. From Cox et al. (in prep.). 

2. Blue supergiants 

Circumstellar structures around BSGs have been predominantely identified as bow 
shocks around runaway stars. Originally discovered with IRAS (e.g. Van Buren & Mc- 
Cray, 1988, ApJ, 329, L93), such structures have also been seen with MSX and WISE 
(Peri et al. 2012). A more general survey of BSGs, i.e. not targeting runaway stars, with 
objects selected from Crowther et al. (2006) and Przybilla et al. (2010), reveals IR ma- 
terial around six of the 45 targets at 22 jam with WISE, also mostly in the form of bow 
shocks (Wachter, in prep). Several examples of bipolar nebulae around BSGs are also 
known (e.g. Sher 25, Smartt et al. 2002; HD 168625, Smith 2007). However, this material 
could have also been ejected during an LBV phase, since LBVs can exhibit BSG spectra, 
and we will therefore concentrate on the bow shocks. 

Runaway stars have large stellar velocities (above 30kms~^) resulting from dynamical 
interactions in (dense) clusters or from a supernova explosion in a binary system. These 
stars can thus travel at supersonic speeds through the local medium giving rise to "bow 
shocks" as their stellar winds interact with the surrounding medium, which has been 
previously ionised by stellar photons from the hot star (Weaver 1977). The occurrence 
of such bow shocks has been shown to depend primarily on the ISM conditions (Huthoff 
& Kaper 2002). For example, even a runaway star may travel at subsonic speeds in the 
tenuous interior of a superbubble, where the sound speed can be as much as lOOkms"^, 
hence no (detectable) bow shock will be produced in that case. The filling factor of ISM 
with '^sound ^ lOkms"^ is 20% and 75% of 0-stars have velocities ^lOkms"^, so the 
expected fraction of 0-stars with bow shocks is ~15%. This is remarkably similar to the 
values derived from IRAS and WISE observations (Noriega-Crespo et al. 1997, Peri et 
al. 2012). 

Once formed, the size, shape and morphology of a bow shock depends on both stel- 
lar (wind kinetic energy and stellar velocity) and interstellar parameters (density and 
temperature). In particular the ratio v^/^wind indicates whether or not instabilities are 
likely to develop (Dgani et al. 1996), and the stand-off distance between the star and 
the apex of the shock is determined from the pressure balance between the stellar wind 
and the ISM (see analytical results by Wilkin 1996 and simulations by e.g. Comeron & 
Kaper 1998, Blondin & Koerwer 1998). Independent estimates of the wind parameters 
can thus be inferred from bow shocks, which serves as a useful check for atmosphere 
models, but the values are sensitive to the ISM properties, which are not always known 
with precision. 

Currently, a small survey with Herschel-PACS of 5 runaways with known bow-shocks is 
ongoing: a Cam, COph, rCMa, Vela X-1 and BD+43°3654 (Cox et al., in preparation). 
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Figure 2. Left: Interferometric image of the photosphere of Betelgeuse obtained by Haubois et 
al. (2009), showing its inhomogeneous surface brightness. Center: VLT/NACO adaptive optics 
tricolor composite image (RGB=KHJ) obtained by Kerveha et al. (2009), showing the emission 
from a compact molecular envelope. Right: VLT/VISIR image at 10.49 /xm of the dust thermal 
emission obtained by Kervella et al. (2011). North is up. East to the left, and the field of view 
is given in the upper right corner of each image. 

For Vela X-1, the peak emission of the dust emission is co-spatial with the most prominent 
Ha arc seen in the supposed direction of space motion (Fig. [T]): it is concluded that the 
outer shock is radiative, but the inner shock is adiabatic, though some Ha emission 
possibly related to (part of) the inner termination shock is also detected. From the 
analysis of its "puffed-up" bow shock (Fig.jl]), the mass-loss rate of BD+43°3654 (04If) 
was found to be lO~^M0yr~^: this is very high (by 2 orders of magnitude) in view of 
current mass-loss rate estimates of such stars, but the exact value strongly depends on 
ISM density, which need to be refined. The dust temperature, ~ 45 K, is compatible with 
heating by stellar photons only, suggesting there is no additional shock-heating of grains. 
The thickness of a bow shock (~ 1 pc) suggests a Mach number close to unity, implying 
a ISM temperature of 10^ - 10"^ K. 



3. Red supergiants 

Circumstellar structures on scales of a few arcseconds or less around RSGs have been 
revealed through interferometric techniques (e.g. Monnier et al. 2004). Stencel et al. 
(1988, 1989) reported the IRAS detection of resolved shells with typical radii of a few 
arcminutes around RSGs for a significant fraction (25%) of their sample. However, higher 
resolution Spitzer images fail to confirm several of these extended structures (Wachter, 
in prep), indicating that a systematic survey is needed to ascertain the occurrence of 
large scale circumstellar shells around RSGs. 

A few (famous) cases have however been studied in depth. One of these is Betelgeuse, 
a cool (3600K), large (700 R©), rather massive (10-15 M©), luminous (> 10^ L©), and 
nearby (150 pc) star. Because of its distance, Betelgeuse can be probed on almost all 
scales, providing a unique panorama of stellar surroundings (Fig. [2|. Space-based and 
interferometric instruments (e.g. HST, lOTA/Ionic and VLTI/Pionier) revealed the pho- 
tosphere, notably the expected non-uniformities due to large convection cells. Adaptive 
optics imaging in the near-IR (e.g. NACO) and (radio or IR) interferometers unveiled 
the properties of the internal, compact molecular envelope (1-10 R*). Precursors of dust 
have been found there, as well as an extended "plume" reaching 6R* and maybe linked to 
a hot spot on the photosphere. High-res imaging (e.g. VLT/VISIR) shows the envelope 
at intermediate scales (10-100 R*), where the dust forms (a possible signature of silicates 
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Figure 3. Continuum subtracted Herschel SPIRE spectrum of VY CMa from 294/im to 
192/im. Multitude of molecular lines have been detected. From Matsuura et al. (in prep.). 

has been found). At these small and intermediate scales, Betegeuse presents a complex 
circumstellar envelope (with knots and filaments) at all wavelengths, which implies an in- 
homogeneous spatial distribution of the material lost by the star. Finally, at the largest 
scale, IR imagers such as Herschel unveil the cool external envelope (100-10000 R*), 
where a bow shock with the ISM is detected (Cox et al. 2012). 

Herschel has also probed the envelope of other red supergiants (Groenewegen et al. 
2011). Turning in particular to the case of VY CMa (Matsuura et al., in prep.), the 
potential of IR spectroscopy is obvious. Herschel-SPIRE reveals a rich spectrum, with 
a dust continuum and hundreds of lines dues to molecules (one third linked to water. 
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Figure 4. Top left: A WISE color composite of 12 /im (blue; blue contours) and 22 /im (red; 
yellow contours) of WSl, discovered and initially characterized by Gvaramadze et al. (2012). The 
contours for each band help illuminate the morphology of the nebular material, which has an 
overall SE-NW elongation, reminiscent of bipolar structure. Bottom left: Optical photometric 
monitoring since discovery show both the R and I light curves have brightened by about 1 
magnitude over the last year. Arrows indicate times when same night spectroscopy were secured. 
Right panel: Optical spectroscopic monitoring indicates evolution to cooler temperature with 
near disappearance of the He I 5876 A and 6678 A and changes in the Hq, line profile. Figures 
from Stringfellow et al. (in preparation). 



others to CO, CS, SiO,...), which constrain the envelope's properties. For example, the 
isotopic ratio ^^C/^^C is found to be 6.5, in agreement with observations of other RSGs 
but at odds with theoretical predictions which are four times higher at least. Very strong 
emission of submm molecular lines can be explained if a temperature gradient is present 
in the envelope, e.g. because of dust formation at a certain radius. 
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Figure 5. PACS spectrum (central spaxel) of WRAY 15-751 nebula, showing the lines from 
the ionized and neutral regions (from Vamvatira-Nakou et al., submitted). 

4. Luminous Blue Variables 

Because of their spectacular eruptions, LBVs are the most well-known cases of massive 
stars with ejecta. It is not yet certain, however, at what stage (BSG? after a RSG phase 
or not?) this material is ejected, and how (multiple events?). LBVs are rare: in the list of 
Clark et al. (2005), there are only 12 confirmed and 23 candidate LBVs. IR has played 
a key role in recent years. The search, through surveys like MIPSGAL, of round-shaped 
nebulae with luminous central stars resulted in the discovery of many new nebulae: 62 
shells in Wachter et al. (2010), 115 shells and bipolar nebulae in Gvaramadze et al. 
(2010), 416 structures in Mizuno et al. (2010). Many of these nebulae are preferentially 
detected with Spitzer 24/im band, indicating relatively cold material. 

Identifying shell-like structures is only the first step. To ascertain a cLBV status, 
the central object needs to be studied spectroscopically. This was done for many of 
these new detections (c.f., Gvaramadze et al. 2010; Wachter et al. 2010; Stringfellow et 
al. 2012a, b, and in preparation). The classification does not rely on the presence of a 
particular line, but rather on the morphological resemblance of the spectra to spectra of 
known LBVs - while not 100% perfect (some peculiar O and WR stars display similar 
features), this method has the advantage of being simple and rather robust. A more 
definitive answer can be provided through photometric and/or spectroscopic monitorings. 
Indeed, as their name indicate, LBVs should be variable. Near-simultaneous photometric 
and spectroscopic monitoring in the optical (and IR) of about a dozen newly identified 
candidate LBVs has revealed that WSl (discovered by Gvaramadze et al. 2012) is indeed 
a bona fide LBV, presently displaying what appears to be S Dor type variability as shown 
in Fig. [4] (Stringfellow et al. 2012, in preparation). 

IR is also useful in revealing details of particular objects. For example, a Herschel 
survey of LBVs undertaken at Liege yielded as first result a characterization of the 
surroundings of WRAY 15-751 (Vamvatira-Nakou et al., submitted). IR photometry can 
only be explained if the star evolves at constant luminosity and dust grains are Fe- 
rich. Images also revealed the presence of a second shell, about 4 times larger than 
the previously known one, which most probably results from an older eruptive event. 
Considering both structures, there is about 0.075 of dust in the system. Ionized gas 
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Figure 6. Examples of the three stages in WR nebula morphologies: from left to right - bubble 
(WR16), clumpy phase (WR8), and mixed phase (WR35b). From Toala et al. (in prep.). 

is responsible for several forbidden lines observed in the Herschel-PACS range (Fig. [5|, 
which allow a N/0 abundance of about 7 times solar and a mass of ionized gas of 1-2 
M0 (20 times that of dust), to be derived. 

Dust can be well studied in the IR, so this range may provide clues on where dust come 
from in galaxies. Two examples of such feedback were presented in the session: rj Car and 
SN1987A. The latter was observed with Herschel at 100-500/im wavelengths, and 0.4-0.7 
M0 of dust was detected - mostly silicates and amorphous carbon (Matsuura et al. 2011). 
It is thought that this dust come from the explosion, but the role played by previous 
mass-loss episodes, in particular the LBV phase, is not yet clear. For example, about 0.12 
M0 of dust was detected, thanks to 30/im MiniTAO observations, in the famous LBV 
T] Car. Up to 80% of that dust belonged to the torus, hence may not be related to the 
big 1843 event. 

5. Wolf-Rayet stars 

Only a few percentage (4-6%) of Wolf-Rayet stars displays surrounding nebulosities in 
the WISE survey, and most are found around WN stars (Wachter, in prep). 

The morphological classification scheme of WR nebulae proposed by Chu (1981) has 
been revised in this meeting by Guerrero et al. (Toala et al, in prep.) using WISE IR 
images and SDSS or Super Cosmos sky survey Ha images for 35 nebulae associated with 
WRs. Two WISE bands were particularly used: the one at 12/im, which encompasses 
PAH lines and lines of low excitation ions, and that at 22/im, to which thermal emission 
from dust and lines of He I as well as high excitation ions contribute. Three phases are 
defined. In the first one, WR nebulae appear as complete shells or bubbles. It corresponds 
to the star just entering the WR stage, when its powerful wind sweeps up the previous 
slow and dense winds (from e.g. LBV or RSG stages). The second phase is the clumpy 
phase. At that point, the nebulae display knots of gas and dust connected by partial 
shells and arcs. It corresponds to an age of a few 10^ yr, when instabilities break down 
the swept-up shell. The stellar motion through the ISM has an impact on the morphology, 
for example one-sided arc may be sometimes seen. Finally, the mixed nebular phase ends 
the cycle, with no definite morphology nor always a 1-to-l correspondence between optical 
and IR images. It corresponds to the last stage, when the circumstellar nebula begins to 
dissolve into the ISM. 
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Figure 7. Visibility of HD 316285 as a function of wavelength for the three different telescope 
baselines measured with VLTI/ AMBER. Notice the drop in visibility within the Bracket gamma 
line, indicating that the line is more extended than the neighboring K-band continuum. From 
Groh et al. (in prep.). 



6. Studying the close environment of massive stars 

The close environment of massive stars is the "missing link" between the star itself 
and the large circumstellar features. It plays a key role in understanding the mass-loss, 
but it is also difficult to probe directly. 

Emission lines arising in the wind and circumstellar material are a classical way to 
study this region, as well as near-IR excess linked to disk-like features. In this context. 
Be and B[e] stars are targets of choice, and surprises are frequent. For example, Graus et 
al. (2012) found three new sgB[e] in the SMC: they display typical spectra, with forbidden 
lines, but the line strengths as well as the IR excess appear reduced compared to usual 
objects of this class. It suggests that either the disks have less material or less dust than 
usual, or maybe that these stars are transitional objects. Another case intriguingly shows 
the opposite situation: CD— 49° 3441 displays forbidden lines and appreciable IR excess, 
but is a main-sequence Be star away from any star- forming region (Lee and Chen 2009). 
A possibility may be that this star is in fact a weak B[e], rather than a classical Be. 

The environment close to the star can also be studied, directly, by means of interfer- 
ometry, which is usually performed at long wavelengths. Most optical/IR interferometric 
measurements rely on the measurements of "visibilities" , which are directly linked to the 
size of the object. Recently, several massive stars, including nine LBVs, were observed 
with the VLTI (Groh et al., in prep.). Amongst these, HD316285 (Fig.!?)): the recorded 
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MIRC Observations of Rapid Rotators 



Regulus Rasalhague Altair Alderamin Bet Cas 



Monnier et al. 2007; Zhao et al. 2009; Che etal. 2011 



Figure 8. Images of rapid rotators derived from MIRC interferometric measurements. Each star 
shows a bright pole and dark equator, which is caused by gravity darkening effect. Modehng of 
Regulus gives gravity darkening coefficient /3 = 0.19, rather than 0.25 as von Zeipel predicted 
in 1924. 

visibilities implied a size of 0.002" for the source of continuum radiation, and 0.004" for 
the source of the Br7 line. A CMFGEN fit to the spectrum yields a stellar model with 
which one can estimate the wind size in IR, and it agrees well with VLTI observations. 
The asymetric shape of the measured differential phases (red vs blue side) favors a prolate 
shape for the rotating star wind, but it could also be explained by clumps or binarity. 
Since the latter imply in time variability, a monitoring will be needed to ascertain the 
exact nature of the asymmetry. 

While visibilities provide valuable data, "real" images are always more impressive. 
Interferometric instruments such as Michigan Infrared Combiner (MIRC) and PIONIER 
are beginning to provide such data. MIRC was the first to image Altair (Monnier et al. 
2007) and several other rapidly rotating stars as shown in Fig. [sj It has also imaged 
circumstellar disks and multi-object systems. For example, the disk contribution of 5 Sco 
was shown to remain stable during the periastron in 2011 (Che et al. 2012), and the 
mass-exchange in the /3Lyr system can be clearly imaged (Zhao et al. 2008), as well as 
the 3 components of Algol (Baron et al. 2012) or the disk of the eclipsing companion of 
e Aur (Kloppenborg et al., 2010). 



7. Conclusion 

This session has demonstrated the usefulness of studies in the IR in studying the 
environment of massive stars. Recent advances in this domain are notably provided by 
surveys, as they enable discovery of new objects to study, thereby improving the census 
of nebular features associated with hot stars. Furthermore, IR diagnostics unveil the 
properties of these neighbouring nebulosities: morphology, temperature, composition, 
density are the necessary keys paving the way of a better understanding of the mass-loss 
in massive stars. 
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